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We present a theoretical study of three-photon electromagnetically induced absorption for a 
ladder-type three-level atomic system. A probe beam was tuned to the lower line and two counter- 
propagating, linearly polarized coupling beams were tuned to the upper line. The system can be 
modeled with a three- (or five-) level scheme when the polarization directions of the coupling beams 
are parallel (or perpendicular). By calculating the absorption coefficients analytically for the two 
schemes, we found that the corresponding absorption coefficients were identical except for differ¬ 
ent transition strengths, and that the primitive scheme embedded in those schemes was a simple 
four-level scheme. 

PACS numbers: 42.50.Gy, 32.80.Qk, 32.80.Wr 


I. INTRODUCTION 

The atomic coherence generated between dipole- 
forbidden states due to two laser beams plays an im¬ 
portant role in the abrupt enhancement of transmission 
or absorption. The effects resulting from this coher¬ 
ence are called electromagnetically induced transparency 
(EIT) and absorption (EIA) , and have drawn 
much interest for their potential application to atomic 
and laser spectroscopies. Particular applications include 
light storage 3, high-resolution spectroscopy Q , non¬ 
linear optics 8|, and quantum information [9j. As in 
traditional EIA, enhanced absorption is observed when a 
third laser beam is introduced in two-beam spectroscopy 
setup, which displays EIT. In this spectroscopic setup, 
three-photon atomic coherence results in enhanced ab¬ 
sorption, rather than the enhanced transmission that 
aries from two-photon coherence. This enhanced absorp¬ 
tion is known to result from constructive interference 
©[HI- It should be noted that destructive interfer¬ 
ence with three laser beams can be observed for different 
transition geometries 

Many studies on three-photon spectroscopy are re¬ 
ported in the literature. A narrow resonance was ob¬ 
served in a four-level atomic system [T^. Chanu et al. 
reported the conversion from EIT to EIA by forming an 
N-type scheme [l^ . which was subsequently explained 
theoretically by analyzing the existing subsystems [l5l| . 
Ben-Aroya also reported lar ge-c ontrast absorption reso¬ 
nances in an N-type scheme |16l | . An optical clock based 
on three-photon resonance was also reported m. En¬ 
hanced absorption with a standing-wave coupling field 
for a A-type atomic system was observed [l^, [l^ • Zhu et 
al. examined dichromatic EIT and four-wave mixing in 
A-type atomic systems [13, [HI . 

Recently, Moon et al. reported three-photon electro- 
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magnetically induced absorption (TPEIA) in a ladder- 
type atomic system where a probe beam was tuned to 
the lower transition line and two counter-propagating 
coupling beams were tuned to the upper transition line 
[23 |. Calculations within the simple four-level model 
showed that this absorption phenomenon was due to 
three-photon coherence. They also studied the relation¬ 
ship between two- and three-photon coherence [13| • As in 
previous studies [H, |13| , where the two coupling beams 
were linearly polarized in perpendicular directions, three- 
photon coherence was the main cause of the enhanced 
absorption. However, when the two coupling beams 
were polarized in the same direction, it was not obvious 
whether or not three-photon coherence was established. 
Enhanced absorption may also have been observed in this 
case. It should be also noted that TPEIA is different 
from the conventional EIA where the absorption results 
from the transfer of the Zeeman coherence of the ex¬ 
cited state The main cause of the absorption in 

TPEIA is the three-photon coherence generated between 
the ground and intermediate states. 


We here present a theoretical and analytical study of 
TPEIA for two schemes, where the coupling-beam polar¬ 
izations are either parallel or perpendicular to one an¬ 
other. To study TPEIA analytically, these two schemes 
were modeled using three-level and five-level schemes. 
After deriving the analytical form of the absorption co¬ 
efficients in the limit of weak probe intensity, we studied 
the similarities and differences between the two schemes. 
These are the new findings as compared to the previous 
reports [13, HI] as well as accurate analytical solutions 
of the absorption coefficients. This paper is structured 
as follows. Section |H] briefly describes the motivation 
and experimental findings of our work. The theoretical 
method of the calculation is presented in Sec. mu Sec¬ 
tion m discusses the absorption spectra averaged over 
the Doppler-broadened atomic velocity distribution. Fi¬ 
nally, we summarize the results. 
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II. MOTIVATION 
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FIG. 1: (Color online) (a) Simplified schematic diagram for 
the TPEIA experiment, (b) Energy-level diagram for the 
5 <S'i/2 — 5 P 3/2 ~ 5 P 5/2 transition in ®^Rb atoms. The two- 
coupling beams are linearly polarized in (c) parallel or (d) 
perpendicular direction. 


The simplified schematic diagram and the energy-level 
diagram for the experiment considered are shown in Fig. 
1. The probe beam is tuned to the 5 S'i/ 2 (F = 2) — 
5P3/2iF' = 3) transition, and two counter-propagating 
coupling beams are tuned to the hPz/ 2 {F' = 3) — 
5D^/2{F" = 4) transition. As shown in Fig. 1(a), all 
the laser beams are linearly polarized, but two polariza¬ 
tion configurations are considered. In the first (paral¬ 
lel) scheme, the polarization vectors of the probe, the 
first, and second coupling beams are in the x, y, and 
y directions, respectively. In the second (perpendicular) 
scheme, the polarization vector of the second coupling 
beam is switched to the x direction. In Fig. 1(a), the 
double-sided arrows and the dots denote the direction of 
the electric field of the laser beams. Thus, the polariza¬ 
tions of the coupling beams are parallel (perpendicular) 
for the first (second) configuration. When the polariza¬ 
tion direction of the probe beam is chosen as the quan¬ 
tization axis, the energy-level diagrams for the parallel 
and perpendicular polarization schemes are those shown 
in Figs. 1(c) and 1(d), respectively. The solid red, solid 
green, and dotted green lines represent the probe, the 
first, and the second coupling beams, respectively. As 
can be seen in these figures, the transition schemes for the 
two configurations are completely different. In the par¬ 
allel configuration, the two coupling beams are tuned to 
the same transition lines, whereas those beams are tuned 
to the different transition lines in the perpendicular con¬ 


figuration. The analytical calculation of the absorption 
coefficients for the basic units of these two schemes is 
discussed in the next section. 
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FIG. 2: (Color online) (i) An experimental EIT spectrum. 
Typical experimental TPEIA spectra for (ii) the parallel and 
(iii) the perpendicular polarization configurations. 


Figure 2 shows typical experimental results for EIT 
and TPEIA. We refer to previous papers for the details 
of the experimental setup [H, The powers of the 
probe and of the first and second coupling beams were 
8 fj,W, 20 mW, and 20 mW, respectively. The diameter 
of the laser beams was 2 mm. In Fig. 2, we present two 
TPEIA spectra for the parallel and perpendicular polar¬ 
ization configurations, and an EIT spectrum obtained in 
the absence of the second coupling beam. In Fig. 2, the 
frequency of the probe beam was tuned to the resonant 
transition 5Si/2{F = 2) — ^Pz/ 2 {F' = 3), whereas the 
frequencies of the coupling beams were scanned around 
the transition hP^/ 2 {F' = 3) — ^D^/ 2 {F" = 2,3,4). As 
discussed previously [H, , the EIT signal transformed 

into the absorption signal when the second coupling beam 
was introduced. As can be seen in Fig. 2 for the two po¬ 
larization configurations, the two TPEIA spectra for the 
5 S'i/ 2 (F = 2 )- 5 P 3 / 2 (F' = 3 )- 5 £i 5 / 2 (F" = 4) transition 
are very similar. As mentioned above, the interaction en¬ 
ergy level diagrams for these two schemes are very differ¬ 
ent. However, we observed very similar TPEIA spectra 
for the two schemes, as can be seen in Fig. 2. This is the 
motivation of our work. 


III. THEORY 

To study the similarity and difference between the re¬ 
sults of TPEIA for the two polarization configurations 
analytically, we used the two simplified energy-level dia¬ 
grams shown in Fig. 3. The transitions (Rabi frequency) 
of the probe, and of the first and second coupling beams 
are denoted by solid red (Hi), solid green (H 2 ), and 
dotted green lines (H 3 ), respectively. Figures 3(a) and 
3(b) show the basic units of the diagrams corresponding 
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FIG. 3: (Color online) Simple energy-level diagram for (a) three-level and (b) five-level schemes, (c) The energy-level diagram 
with multiple alternative interactions for the upper transition line, equivalent to the three-level scheme, (d) Four-level scheme. 


to the parallel and perpendicular-polarization configura¬ 
tions, respectively. In Fig. 3(b), (f^a) denotes the 

Rabi frequency of the first (second) coupling beam for 
the other transition line. Figure 3(c) shows the energy- 
level diagram with multiple alternative interactions for 
the upper transition line. This diagram was proved to 
be equivalent to the schemes in Figs. 3(a) and (b) and 
will be discussed below in detail. Figure 4(d) shows the 
four-level scheme, the simplest primitive scheme for the 
two schemes in Figs. 3(a) and (b). 

We calculated analytically the absorption coefficients 
of the probe beam in the weak-probe-intensity limit, for 
an arbitrary intensity of the coupling beams. We first 
considered the three-level scheme shown in Fig. 3(a). 
The decay rates of the states |1), |2), and |3) are Fi, F 2 , 
and r 3 (= 0 ), and the decay rates of the optical coher- 
ences between the states |/x) and \v) are 7 ^ 1 , = 
with < V. The wavelengths, wave vectors, and detun¬ 
ings of the probe (coupling) beam are Ai, fci, and Sp (A 2 , 
k 2 , and 5c), respectively. The effective detunings (experi¬ 
enced by an atom moving at velocity v) of the probe and 
of the first and second coupling beams are i5i = 5p — kiv, 
S 2 = 5c + k 2 V, and 63 = 5c — k 2 V, respectively. 

We considered the solutions up to the first order in the 
probe Rabi frequency. In the energy level diagram shown 
in Fig. 3(a), the two coupling beams were tuned to the 
transition between |1) and |2). Therefore, the frequency 
mixing produced a variety of oscillation frequencies for 
the optical coherences and populations. Because of the 
weak probe intensity, the optical coherence pi 2 can be 
neglected, and thus we need consider only the optical co¬ 
herences P 23 and P 13 . In addition, all the stationary and 
oscillating terms in the populations vanish except for the 
stationary term of P 33 , which is unity in the limit of weak 
probe intensity. If the probe beam is not quite weak, 
the coherent population oscillation (CPO) in the popu¬ 
lations of the ground, intermediate, and excited states 
should be considered . In this case we would need 

a numerical calculation instead of the analytical calcula¬ 
tion. We refer to previous publications for the method of 
calculation of the oscillation frequencies [26l-l^. When 


Wphoton interactions are taken into account, the optical 
P23 and pi 3 are explicitly given by 

P 23 = ui+ U 3 e*'^'^‘ -b V 3 e~''^'^* + -b 

-I-b -b (la) 

P13 = U2 + V2e~"^'^* + -b H- 

+ (lb) 

where 5^ = 53 — 52 = —2k2V and N is assumed to be an 
odd integer, u and v in Eqs. (Hal) and m are Fourier 
components, which should be obtained by solving the 
density matrix equations, and the subscripts of u and v 
represent the number of photon interactions. Thus, the 
Fourier components of P 23 (pia) result from an odd (even) 
number of photon interactions. From the density matrix 
equations, we obtained the differential equations for the 
components of P 23 and P 13 as follows. The lowest-order 
equation for iii is given by 

2 2 2 

ill = iAiUi - -Hi - -Q. 2 U 2 - -^3V2- (2) 

We can see that the equations for the u- and u-series can 
be constructed separately. The equations for the u-series 
are given by 

U 2 = f A 2 U 2 — ipiUi — ip2U3, 

U3 = f A3U3 - ip2U2 - ip3UA, 

iln = '^Pn—l'^n—1 ‘^Pn'^n+1, (3) 


and those for the u-series by 

V2 = iA2V2 - iqiUi - iq2V3, 

V 3 = iA'^Vs - iq 2 V 2 - iq3VA, 

Vn = iA'^Vn - iq-n-lVn-l “ iqnVn+1, (4) 
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where vi = ui, and 


final result for ui: 


Pn{qn) 


% (%) : for ri= 1,3,5, ••• 
%(^) , for n = 2,4,6,--- . 


( 5 ) 



The effective detunings are given by 

A2„ = di + n(52 - (n - 1)^3 + *713, 
A2„+i = (5i + n<52 - + *723, 

A 2 n = (5i + nda - (n - 1)^2 + *713, 
A2 „+i = di + nda - n(52 + *723- 


Equations ©“(H]) can be described graphically as in Fig. 
3(c). As we discuss below, the results for the schemes 
in Figs. 3(a) and (b) are equivalent to the energy-level 
diagram in Fig. 3(c). The right and left branches denote 
the equations for the u- and u-series, respectively. Since 
these equations are equivalent under exchange of 112 ( 1 ^ 2 ) 
and r 23 (< 53 ), we need only solve the solutions for either 
the u- or u-series. Here, we solve the former. 

We assumed the following relation: 


Pn-l 

Un — . Ufi — l’ 

■^n 


( 6 ) 


If we insert Eq. m into Eq. ([S]) by replacing n with 
n -I- 1, i.e., Mn-i-i = obtain 


Un 


Pn-l 


A, 



^n—1 • 


( 7 ) 


Comparison of Eqs. and © yields the following re¬ 
cursion relation: 


A„ = A 


n 


P 


2 

n 


An-t-l 


( 8 ) 


Because the maximum value of n is A*, Ajv = A^r and 
An = 0 when n > N. 

Using Eqs. ®, dH), and (|S]) , we can express U 2 in terms 
of ui as 


-1 



Similarly for U 2 , 


-1 



Inserting Eqs. ® and (nni into Eq. ®, we obtain the 
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Hi/4 
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A 3 —- 




A 4 . 
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A'-- 


ni/4 


ai/4. 


Ai-- 


( 11 ) 


where ui for the parallel polarization configuration is de¬ 
noted U||. We notice that the terms associated with A 3 
and A 3 are responsible for TPEIA and result from the 
three-photon coherences U 3 and V 3 in Eq. (ITal) . respec¬ 
tively. As shown in the next section, when thermal aver¬ 
aging is considered, the contribution of U 3 is much larger 
than that of V 3 owing to the Doppler-free characteristics. 
It should be noted also that the term U 3 in Eq. (fTal) . 
responsible for the four-wave mixing signal, is given by 



( 12 ) 


We performed the calculation for the five-level model in 
Fig. 3(b), the basic unit for the energy-level diagram for 
the perpendicular polarization configuration. As in the 
three-level scheme, we considered the solutions up to the 
first order in the probe Rabi frequency. Thus, we need 
only consider the elements P35, /945, P15, P25, and P55(— 1 ) 
while neglecting all the other elements. From the calcu¬ 
lation of the oscillation frequencies, these density-matrix 
elements can be expanded as follows: 


P35 

= ui + -1- Use 




(13a) 

P45 

= U3e*'*‘^* -1- 



-1- -1 -, 

(13b) 

P 25 

= U2 -1- 



-h H -, 

(13c) 

Pl5 

= -f 



-1- -. 

(13d) 


Comparing Eqs. (I13a|) - (ll3dl) with Eqs. (flail and (fib)) , it 
is easy to recognize that the terms in P 35 and P 45 are the 
same as those in P 23 in Eq. (Hal) . In addition, the terms 
in P 25 and pis are the same as those in pi 3 in Eq. m. 
Thus, the structures of the density matrix elements of 
the three- and five-level schemes are identical. 

When we solve the density-matrix equations using the 
expanded elements in Eqs. H3al) - H3dl) . we obtain equa¬ 
tions identical to those shown above in the discussion 
of the three-level scheme. The only difference is in the 
definitions of and because of the different Rabi fre¬ 
quencies that are due to the different transition strengths. 
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In this case, Pn and qn are defined as 


Pn{qn') — ^ 



^2 ( t 
2 V 2 j ’ 


2 V 2 / ’ 


for n = 1, 5,9, • • • , 
for n = 2 , 6, 10 , • • • , 

for n = 3, 7,11, • • • , 
for n = 4, 8, 12, • • • . 


Then, the final value of ui for the perpendicular config¬ 
uration is given by 


Hi 



where ui for the perpendicular-polarization configuration 
is denoted u±. Equation m is equivalent to Eq. m 
except for the different values of the Rabi frequencies. It 
should be noted that the three-photon coherence in 
Eq. (I13bl) is responsible for TPEIA. 

As discussed in the next section, the terms of the v- 
series in Eqs. (El) and ED are not Doppler-free, so that 
the most primitive scheme for the parallel and perpen¬ 
dicular schemes is the simple four-level scheme shown in 
Fig. 3(d). The result for this scheme is simply given by 


-1 



where U 4 l is ui for the four-level scheme. 

Finally, the absorption coefficient for each polarization 
configuration averaged over the Doppler broadened ve¬ 
locity distribution is given by 


1 SAf iVatr2 


V^Uinp StT H 


-p pOO 

— / du Im(Mj)(16) 

1 J — OO 


where j G {||j-L, 4L}, Ump is the most probable speed, 
and A^at is the atomic density in the vapor cell. 


IV. RESULTS AND DISCUSSION 

The calculated normalized absorptions for a stationary 
atom in the parallel-polarization configuration with var¬ 
ious values of N are presented in Fig. 4. When N — 1, 
the absorption is the background value in the absence of 
the two coupling beams. We obtained an EIT spectrum 
for N = 2 where the second coupling beam is absent, 
and a TPEIA spectrum for V = 3 in the presence of two 




FIG. 4: (Color online) Calculated normalized absorptions for 
a stationary atom at (a) A = 1,2, and 3 and (b) N = 1,4, 
and 12. 


coexisting coupling beams. In Fig. 4, as N increases, 
the transmittance and absorption vary in an alternating 
fashion [l^. As will be seen later, there is a big differ¬ 
ence between the results for TV = 2 and N = 3, but no 
significant differences for the higher values of N. 

In order to study the effect of thermal averaging, 
we present the absorption coefficient for the four-level 
scheme depicted in Fig. 3(d). We used this scheme 
because this is the most primitive in the exact energy- 
level diagrams of the parallel and perpendicular configu¬ 
rations. In addition, it can simplify the occurrence of ab¬ 
sorption resonances. The four-level scheme corresponds 
to A = 3 where only the u-branch exists. This choice 
will be validated soon. Figure 5 plots the normalized ab¬ 
sorptions for the velocities v = 0, ±5, and ±10 m/s, and 
the thermal-averaged absorption. The resonances occur 
at = 0 and ±(l/ 2 )y^H 2 ± when u = 0 , from a 
dressed-state analysis. As the magnitude of the veloc¬ 
ity increases, the two side resonances shift far away from 
the origin, but the central resonance remains almost sta¬ 
tionary at the origin. Therefore, after averaging over the 
velocity distribution, the central peak survives, while the 
two side resonances cancel out. Although the analysis 
was performed for the four-level scheme, this analysis is 
valid for the other complicated diagrams in Fig. 3. 
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Coupling Detuning (MHz) 


FIG. 5: (Color online) Calculated absorptions for an atom 
with the four-level scheme, moving at velocity u = 0, ±5, or 
±10 m/s. The bottom curve refers to the thermal-averaged 
absorption spectrum. 


Figure 6 shows the series of averaged absorption spec¬ 
tra for the analytical results in Eq. (HU). Figure 6 (a) 
shows the absorption spectra for several values of the 
Rabi frequency of the second coupling beam, where 
rii = 27r X 1 MHz and 172 = 27r x 10 MHz. As can be seen 
in Fig. 6 (a), the EIT spectrum at fla = 0 transforms to 
the TPEIA spectrum as fla increases. Figure 6 (b) shows 
three typical spectra for the parallel-polarization config¬ 
uration where = 27r x 1 MHz, 172 = ^3 = 27r x 10 
MHz, and the probe detuning was set to zero. The black 
curve (A) denotes the exact result of Eqs. (ITTI) and (fT6l) . 
The red curve (B) denotes the spectrum where the term 
of the u-series was neglected. Clearly, the two spectra 
are very similar. This means that the contribution of the 
u-series is much smaller than that of the u-series. The u- 
and u-series are represented graphically as the right and 
left branches in Fig. 3(c). The reason for the weak con¬ 
tribution of the V series is that it is not Doppler-free. In 
Eq. (HH), the two-photon resonance terms for the u- and 
u-series are A 2 = ± (52 ± *713 and A 2 = ± (53 ± 1713 , 

respectively. Because ± (52 = (5p ± (5c ± {k 2 — ki)v and 
<5i ± <53 = (5p ± 5c — (^2 ± ki) V, the term A 2 contributes 
much more than A 2 after averaging over the velocity dis¬ 
tribution. The blue curve (C) in Fig. 6(b) shows the 
result for the four-level scheme depicted in Fig. 3(d). 
This is very similar to the exact calculated spectra. As 
172 and Da decrease, the discrepancy between the four- 
level and exact results decrease. Figure 6(c) shows the 
same spectra as in Fig. 6(b), but with the probe detun¬ 
ing being scanned while the coupling detuning is set to 
zero. The overall spectra in Figs. 6(b) and (c) show very 
similar trends. 

Finally we consider the conditions for detunings of the 
two coupling beams to create TPEIA signal, which were 
set to zero in the paper. Now we assume that the coupling 
detunings are arbitrary, and thus S 2 = Sci + k 2 V and 
^3 = Sc 2 — k 2 V where 5ci and 5c2 are the detunings of the 
first and second coupling beams. Equation (TTSl) can be 



Coupling Detuning (MHz) 



FIG. 6: (Color online) (a) TPEIA spectra for various Da val¬ 
ues. TPEIA spectra obtained via three different calculating 
methods with scanning of (b) coupling detuning and (c) probe 
detuning. 


written as 

_ 4A2A3 - D| _ 

2 4AiA2A3-(A3D2±AiD2)- ^ > 

In Eq. (II7L we have two conditions for the two- 
photon Doppler-free characteristics such as A 2 = 0 and 
A 3 DI ± AiD§ = 0. These conditions result in the re¬ 
quirements 6 = —5cl and 6 ci = — 5 c 2 when ki ~ ^2 and 
D 2 ~ D 3 . Therefore, the frequencies of the first and sec¬ 
ond coupling beams should be located at the symmetric 
positions relative to the upper resonant transition line. 
In our case, because 5ci = Sc 2 = 0, we see the TPEIA 
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signal at the zero probe detuning as shown in Fig. 5. 

V. CONCLUSION 

We have presented a theoretical study of three-photon 
electromagnetically induced absorption in ladder-type 
atomic systems. The primitive energy-level diagrams for 
the parallel and perpendicular coupling-beam polariza¬ 
tion configurations were the three- and five-level schemes, 
respectively. We derived analytically the absorption co¬ 
efficients for these two schemes and found them to be 
equivalent to a ladder system with an upper transition of 
multiple alternative interacting schemes. The absorption 
coefficient consisted of contributions from the two oppo¬ 
site branches, one of which was Doppler-free. It was also 
found that the simplest scheme was the simple four-level 
scheme. 

In contrast to normal EIA in degenerate two-level sys¬ 
tems, where the transfer of coherence or population plays 


a major role 1291. 3011 . TPEIA involves constructive in¬ 
terference [Ifl, [ill three-photon coherence is the 
main cause of TPEIA in both parallel and perpendic¬ 
ular schemes. Although the basic mechanism of TPEIA 
emerged from the analytical calculation of the absorp¬ 
tion spectra, TPEIA for real atoms displayed interesting 
behavior. For example, the polarization dependence was 
different depending on whether or not the transition was 
closed. An accurate calculation and analysis for TPEIA 
for real atoms are currently in progress. 


Acknowledgments 

This research was supported by Basic Science Re¬ 
search Program through the National Research Foun¬ 
dation of Korea(NRF) funded by the Ministry of Sci¬ 
ence, ICT and future Planning(2014RlA2A2A01006654 
and 2015R1A2A1A05001819). 


[1] S. E. Harris, Phys. Today 50, 36 (1997). 

[2] F. Fleischhauer, A. Imamoglu, and J. P. Marangos, Rev. 
Mod. Phys. 77, 633 (2005). 

[3] A. M. Akulshin, S. Barreiro, and A. Lezama, Phys. Rev. 
A 57, 2996 (1998). 

[4] A. Lezama, S. Barreiro, and A. M. Akulshin, Phys. Rev. 
A 59, 4732 (1999). 

[5] A. V. Taichenachev, A. M. Tumaikin, and V. I. Yudin, 
Phys. Rev. A 61, 011802(R) (1999). 

[6] M. D. Lukin, Rev. Mod. Phys. 75, 457 (2003). 

[7] A. Krishna, K. Pandey, A. Wasan, and V. Natarajan, 
Eurphys. Lett. 72, 221 (2005). 

[8] D. A. Braje, V. Balic, S. Goda, G. Y. Yin, and S. E. 
Harris, Phys. Rev. Lett. 93, 183601 (2004). 

[9] K. Hammerer, A. S. Sprensen, and E. S. Polzik, Rev. 
Mod. Phys. 82, 1041 (2010). 

[10] N. Mulchan, D. G. Ducreay, R. Pina, M. Yan, and Y. 
Zhu, J. Opt. Soc. Am. B 17, 820 (2000). 

[11] L. Zhang, F. Zhou, Y. Niu, J. Zhang, and S. Gong, Opt. 
Commun. 284, 5697 (2011). 

[12] M. Yan, E. G. Rickey, and Y. Zhu, Phys. Rev. A 64, 
043807 (2001). 

[13] M. G. Bason, A. K. Mohapatra, K. J. Weatherill, and G. 
S. Adams, J. Phys. B 42, 075503 (2008). 

[14] S. R. Chanu, K. Pandey, and V. Natarajan, Eur. Phys. 
Lett. 98, 44009 (2012). 

[15] K. Pandey, Phys. Rev. A 87, 043838 (2013). 

[16] 1. Ben-Aroya and G. Eisenstein, Opt. Express 19, 9956 

( 2011 ). 

[17] T. Hong, C. Gramer, W. Nagourney, and E. N. Fortson, 


Phys. Rev. Lett. 94, 050801 (2005). 

[18] S. A. Babin, D. V. Ghurkin, E. V. Podivilov, V. V. 
Potapov, and D. A. Shapiro, Phys. Rev. A 67, 043808 
(2003). 

[19] I. H. Bae, H. S. Moon, M. K. Kim, L. Lee, and J. B. Kim, 
Opt. Express 18, 1389 (2010). 

[20] J. Wang, Y. Zhu, K. J. Jiang, and M. S. Zhan, Phys. 
Rev. A 68, 063810 (2003). 

[21] G. Q. Yang, P. Xu, J. Wang, Y. Zhu, and M. S. Zhan, 
Phys. Rev. A 82, 045804 (2010). 

[22] H. S. Moon and T. Jeong, Phys. Rev. A 89, 033822 
(2014). 

[23] Y. S. Lee, H. R. Noh, and H. S. Moon, Opt. Express 23, 
2999 (2015). 

[24] J. Fuchs, G. J. Duffy, W. J. Rowlands, A. Lezama, P. 
Hannaford, and A. M. Akulshin, J. Phys. B 40, 1117 
(2007). 

[25] T. Laupretre, S. Kumar, P. Berger, R. Faoro, R. Ghosh, 
F. Bretenaker, and F. Goldfarb, Phys. Rev. A 85, 
051805(R) (2012). 

[26] H. R. Noh and W. Jhe, Phys. Rev. A 75, 053411 (2007). 

[27] H. R. Noh, S. E. Park, L. Z. Li, J. D. Park, and C. H. 
Cho, Opt. Express 19, 23444 (2011). 

[28] S. E. Park and H. R. Noh, Opt. Express 21, 14066 (2013). 

[29] C. Goren, A. D. Wilson-Gordon, M. Rosenbluh, and H. 
Friedmann, Phys. Rev. A 67, 033807 (2003). 

[30] C. Goren, A. D. Wilson-Gordon, M. Rosenbluh, and H. 
Friedmann, Phys. Rev. A 70, 043814 (2004). 




